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The critical behavior of the Ising chain with long-range ferromagnetic interactions 
decaying with distance r", 1 < a < 2, is investigated using a numerically efficient 
transfer matrix (TM) method. Finite size approximations to the infinite chain are 
considered, in which both the number of spins and the number of interaction con- 
stants can be independently increased. Systems with interactions between spins 
up to 18 sites apart and up to 2500 spins in the chain are considered. We obtain 
data for the critical exponents v associated with the correlation length based on 
the Finite Range Scaling (FRS) hypothesis. FRS expressions require the evaluation 
of derivatives of the thermodynamical properties, which are obtained with the help 
of analytical recurrence expressions obtained within the TM framework. The Van 
den Broeck extrapolation procedure is applied in order to estimate the convergence 
of the exponents. The TM procedure reduces the dimension of the matrices and 
circumvents several numerical matrix operations. 



1 Introduction 



It has been well-known that the thermodynamical properties and critical behavior of phys- 
ical models is affected by the presence long-range interactions. They define completely 
new classes of models, where cooperative effects are enhanced and change the thermody- 
namical properties in comparison to those of the short range interaction models. Due to 
the presence of simultaneous coupling among many degrees of freedom, the long range 
models offer very difficult technical problems that makes it impossible to derive exact ex- 
pressions for the thermodynamical properties within the equilibrium statistical ensemble 
formulation. This is observed already for the most simple long-range Ising chain, which 
has been investigated for many decades. In its most simple version, each spin ai interacts 
with all other spins on the chain mediated by coupling constants Jr = J/r"", where r is 
the distance between the interacting spins measured in integer number of lattice spacings. 

Although no closed form solution for this model is available, there are several rigorous 
results on the existence of distinct thermodynamical phases, what depends on the range 
of values of a. The most important features of these rigorous results are: for a > 2, the 
system shows only a disordered phase, VT [T1I2]; for 1 < a < 2, there is a phase transition 
at finite temperature Tc OEj; critical mean field behavior occurs for 1 < a < 1.5 Pfl |2]; 
for a < 1, only one ordered phase exists, VT. Regarding the evaluation of approximate 
results, both renormalization group schemes [Sj-[H] and numerical calculations of finite size 
systems [Zj-fOj have been used to estimate the thermodynamical properties, the critical 
temperature and the critical exponents when 1 < a < 2. 

More recently, the scaling behavior of this model, specially in the range a < < 1, 
where the energy is not an extensive quantity, has also been addressed. Investigations have 
been motivated by an universal scaling scheme proposed by Tsallis that should 

be valid for both extensive and non-extensive long range models, as those constituted by 
spins |I1|-|IZI rotors [IHllin], and so on. 

In this work we use a numerically efficient transfer matrix (TM) approach [SHj to 
analyse the critical behavior of the system, estimating the critical exponents u associated 
with the correlation length in the range for 1 < a < 2, where the system undergo a phase 
transition. In a previous paper [21] , we have introduced this approach to check the validity 
of Tsallis' scaling conjecture to the Ising long range chain, and to find estimates for the 
critical temperature. We obtained results that show very good accordance with other 
numerical estimates, indicating that the proposed approach is quite reliable. It takes into 
account the long range interactions between spins up to a certain distance g apart in the 
evaluation of the thermodynamical properties of the system. Also, this method allows to 
independently increase the number of spins N in the chain, so that N > g + 1. 

The TM approach used herein makes use of very compact matrices, so that the con- 
figuration energies and Boltzmann weights, that are numerically evaluated, can be stored 
and operated in a very efficient way. For instance, we note that the TM procedure avoids 
the numerical evaluation of the TM eigenvalues. For the evaluation of z/, this framework 
is quite useful, as it leads to analytical expressions for the derivatives of the Boltzmann 
weights that can be similarly stored in very compact matrices. So, the derivatives of 
thermodynamical properties can be directly computed, avoiding the use of numerical 
differentiation [7j. 
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Figure 1: A schematic representation of the construction of the infinite 3-chain [g = 3) 
for c = 0, 1 and 2. 

This work is organized as follows: in Section II we discuss the essential aspects of 
transfer matrix method (TM) used to evaluate the critical exponents. In Section III 
we present the critical temperature, the critical exponents u and f3 associated with the 
correlation length and the magnetization for values of a in the following range: 1 < a < 2. 
In Section IV we discuss the results and compare with some previous results [Zj. Section 
V closes the work with our final remarks and perspectives. 

2 Transfer matrix framework 

The long range Ising chain is described by the Hamiltonian 

oo oo oo 

T-(- = J2 JrOTiai+r - hJ2'^ii (1) 

i=0 r=i+l 1=0 

One can construct this model by starting with an infinite chain but a finite number 
interaction constants. Then, in a series of steps g, all Jg are introduced at once, among all 
pairs of spins that are g sites apart. It is also possible to construct the model by taking a 
finite chain with g distinct coupling constants and N = g + c+1 > g + 1 spins. Thus, for a 
generic value g, ii c new spins are included in the chain, g x c new interactions, mediated 
by Jr, r = 1,2, ...g, are added, linking every new spin to all spins up to g sites apart that 
have been introduced at the previous generations. In Figure la, lb and Ic, we show, for 
example, for fixed g = 3, the interactions for the first 3 values of c. If g oo, oo, 
so that this procedures leads to the same the actual system. However, for finite values of 
g, the first construction procedure can be recovered by letting c — > oo. 

In this work, we consider the second way of construction quoted above, which is 
schematically represented in Figure 1, for which the finite chain has N = g + c + 1 spins 
that interact with all spins g sites apart. This way, it is possible to obtain a partition 
function based on 2 x 2 TM Mg^c, instead of 2^ x 2^ TM Mg introduced in reference [7]. 
Besides working with smaller TM's, this approach has the advantage of avoiding the use 
of numerical diagonalization procedures to evaluate the largest eigenvalues. According to 
this TM scheme, which has been described with enough details in j2I], the TM 2x2 TM 
Mg^c are given by: 
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M, 



g,c 



{l[Pk){Qg-Pgr-Lg = .Rg,,Lg. 

k=l 



where Lg is a 2^^'^ x 2 matrix whose elements {Lg)ij = 1, for i + j even, and {Lg)ij 
for i + j odd; Pk are recursively expressed by: 



(2) 
= 0, 



{Pk)i,j 



(p ^ ^(-1)'"' 



for i < 2'=-\ 
and j < 2^ 

■fcj2fe-j+i,2fe+i-j+i for 2 < z < 2 

and 2^' < j < 2^=+^ 
otherwise 



(3) 



where = exp{Jk/T)] Qg is a 2^+^ x 2^ matrix defined by 



1 for i = j or i = j + 2^ 
otherwise 



(4) 



Therefore the free energy per spin, fg c = —T\n{Zg c)/N follows from the partition 
function 

Zg,, = 2A+, = 2((M,,,)i,i + (M,,,)i,2) = T.iRs,ckr (5) 



The correlation function between the first and the r-th spins along the chain, restricted 
to the case c = 0, defined by a (^-dependent correlation function Cg{r;T) = {(Ticrr)g,r = 
1, ...,g, is given by 

1 



Cg{r;T)^-Y: 



9 i,j 



9 

k=2 



P,{X{Pk)Lg,r 



where 



El 



ai 6i 
-bi -ai 



29- 



(6) 



(7) 



with L[x] = largest integer in x. 

The definition of Cg{r;T) can be extended to Cg^i.i + r;T), where c > 0, r > (7. If 
we use (0) for r = g, we have Cg{g] T) = A~/A+. Note further that the correlation length 
for a chain composed of patches described by the matrix Mg c=o is given by 



ln(A+/A, 



HCgig;T,,g)y 



Based on the above expressions for thermodynamical properties, we obtain analytical 
expressions for A+ and A^ given by: 



29- 



A 



± 



5:[a(j,l)±a(2^ + l-j,l)]7,^(j)], 



(9) 



where 
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7g^+i = a(2j-l,2)7±(2j-l) + a(2j,2)7±(2j) j = 1, 2, . . . , 2^- 
y±+i = ±a(2j-l,2)7± 



7± = ±a(2j-l,2)7±(2^' + 2-2j)±a(2j,2)7±(2^' + l-2j), j = 2^'-^ + 1 . . . 2^'-^ 



(10) 

with 7±(0) = 0. 

Since the eigenvalues are given by 0, we can easily derive expressions for their 
derivatives with respect to T and from which the thermodynamic functions, expressed 
in terms of the derivatives of the free energy, can be easily evaluated without resorting to 
numerical differentiation. 

To obtain the scaling properties of ^ we will make use of the FRS framework [7] . This 
scheme proposes a scaling hypothesis which is formally similar to the well known finite 
size scaling, which compares the behaviors of finite size systems with different number of 
components. In FRS, one assumes similar relations among systems with distinct number 
of coupling constants, hence of different interaction ranges. Starting from the assumption 
that, close to the critical temperature T, any thermodynamical function y{t) of the reduced 
temperature t = {T — T/T) is described by a power law 

y{t) = Aot-P, (11) 

and that a finite number g of coupling constants modifies the actual criticality by a 
correction factor /, it is proposed that 

yg{t) = yoo{t)fig/U)- (12) 
With (fTTj) and (fT^ it is possible to show that, close to T, the condition 

Ut) _ 9 



holds. Or, equivalently, t and t' are related by 



(13) 



t'-l^]it. (14) 



The critical temperature Tg at order g is obtained by the condition 

^aiig) _ 9 



(15) 



^g+i{tg) g+V 

and it is expected that the series of values tg converges to the actual t = in the limit 
g ^ oo. _ 

Linearizing and expanding around Tg, leads to 

e,+ l(t) ^g+^{tg) + id^g+^itg)/dt)t 



Ut') Ut9) + idUt9)/dt)t' 

Combining (fT^ with (fTHjl . we obtain estimates of the critical exponent Ug as 



(16) 



] / d(g(tg)/dt X _ , 
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Figure 2: Dependence of c on c, for a = 1.8 and a fixed value of g, for example, g = 10. 

3 Results 

In reference j,21J, we have obtained a series of estimates T for the critical temperature T, 
based on the observed behavior for Cg{r;T), which corresponds exactly to the condition 
expressed by (fT3j). In such case, we have considered N = g + 1, i.e., c = 0. The actual 
value of T has been evaluated by an extrapolation procedure, with very good accuracy, 
although the individual values Tg are not so close to T. 

To address the question of the critical behavior near critical temperature, we need to 
be in a much closer neighborhood of T. So we found it becomes necessary to consider a 
larger number of spins in the chain, what requires larger values of g and c. Increasing g 
leads to exponential growth of storage capacity and computing time. We found, however, 
that taking a better approximations to the infinite chain, what amounts to increasing c 
for a fixed value of g, makes it possible to reach a closer neighborhood of T. 

To this purpose, we evaluate a series of critical temperatures Tg c depending on both 
g and c, from which we obtain first Tg = limc^oo^g,c and T = lim^^oo^g- Each Tg c is 
obtained by imposing the condition (fTBj). where each ^g actually indicates ^g^c, which is 
evaluated for corresponding values of g and c with the help of In Figure 2, we show 
how Tg^c depends one, for a fixed value of g when a = 1.8; the behavior does not change 
for other values of g. 

In our numerical evaluation, Tg has been obtained with a given accuracy, as we increase 
the value of c until | Tg c+i — Tg ^ \ < e, where the highest accuracy corresponds to e ~ 10~^. 
For such situation, c can be as large as 2500, depending on the values of g and a. The 
largest value of = 18 has been obtained for a = 1.6. 

In Tables 1, 2, 3 and 4, we show the values of Tg and Ug, for distinct values of g and 
1 < a < 2: a = 1.2, 1.4, 1.6 and 1.8, for a fixed value of accuracy e = 10~^. In order 
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a = 1.2 




f — 10-^ 






f — 10-3 

c — X \J 




y 


r(f) 


T 






T 


'^9,c 




684 


5 5419339 


2 4783289 


208 


5 3921 944 


9 433681 1 




1058 


6 6798492 


2 7505706 


314 


fi 4445671 


9 6616931 


7 


1 3Qfi 


7 3835491 


2 9643698 


41 


7 0798148 


9 831 5859 


9 


1706 


7 8632795 


3 1338069 


492 


7 4796335 


2 9552605 


11 


1990 


8.2103090 


3.2697569 


566 


7.7599961 


3.0455192 


13 


2260 


8.4726154 


3.3805739 


632 


7.9590277 


3.1099467 


15 


2500 


8.6753271 


3.4710018 


692 


8.1041143 


3.1550949 


VBS extrapolation 




9.482310 


3.803537 




8.910613 


3.288895 



Table 1: Values of Tg and Vg, for odd values of g and a = 1.2. Data are displayed for 
accuracy e = 10~^ and e = 10"^ that was required for the limiting value of Tg. The 
corresponding smallest value of c at which the accuracy was reached is also indicated. 



a = 1.4 




e = 10-^ 
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C(6) 


T 




3 


492 


3.8408632 


2.1731057 


5 


694 


4.3296422 


2.2598101 


7 


856 


4.5894104 


2.31620273 


9 


994 


4.7474461 


2.35392944 


11 


1112 


4.8515617 


2.37937145 


13 


1216 


4.9241076 


2.39645998 


15 


1310 


4.9767888 


2.40768914 


VBS extrapolation 




5.022703 


2.437089 



Table 2: Values of Tg and Ug, for odd values of g and a = 1.4. Data are displayed for 
accuracy e = 10—^ that was required for the limiting value of Tg and the corresponding 
smallest value of c is also indicated. 

to analyse the effect of the accuracy e, for values of a, inside and outside the mean field 
region, a = 1.2 (Table 1) and 1.8 (Table 4) respectively, we include, in these tables, the 
data for another value of e, e = lO"^. It is easy to see that the value of c is very sensible 
to the accuracy of the corresponding values of c for distinct values of e; in Figure 2 we 
plot the critical temperature for increasing values of c, corresponding to, respectively, 
decreasing values of e. 



We observe that the estimates Tg c, for any value of a, are much larger than those 
obtained previously for c = that we have obtained in[23- We observe that, despite the 
estimates c are much larger than c=o; some of the extrapolated values are smaller 
than those that were predicted from the series c=o if the value of e is not small enough. 
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rv — 1 




f — 10-4 




n 

y 




T 






389 


2 81 78502 


9 0697087 




51 n 


3 0448535 


9 0947359 


7 


fins 


3 1 555507 


9 1059548 


q 


fi88 


3 91 90649 


9 1 091 374 




75fi 


3 9591 533 


9 1095755 


13 


814 


3 2859512 


2 1080216 


15 


866 


3.3047686 


2.1053158 


17 


910 


3.3180850 




18 


932 


3.3234981 




VBS extrapolation 




3.328188 


2.109820 



Table 3: Values of Tg and z/^,, for odd values of g and a = 1.4. Data are displayed for 
accuracy e = 10"^ that was required for the limiting value of Tg and the corresponding 
smallest value of c is also indicated. 



a = 1.8 




e = 10-4 






e = 10-3 
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c{e) 


T 




c(e) 


T 


^9,c 


3 


314 


2.1323856 


2.1231430 


94 


2.0633542 


2.1560530 


5 


408 


2.2243788 


2.1639020 


122 


2.1346603 


2.2165373 


7 


482 


2.2642364 


2.1873097 


142 


2.1575740 


2.2647878 


9 


540 


2.2847804 


2.2036550 


158 


2.1642335 


2.3109579 


11 


592 


2.2967828 


2.2163352 


172 


2.1644328 


2.3583671 


13 


636 


2.3039057 


2.2273095 


186 


2.1629377 


2.4063897 


15 


676 


2.3083424 


2.2372586 


198 


2.1588773 


2.4598516 


VBS extrapolation 




2.315334 


2.246300 




2.167883 


2.910225 



Table 4: Values of Tg and Ug, for odd values of g and a = 1.8. Data are displayed for 
accuracy e = 10~4 and e = 10"^ that was required for the limiting value of Tg and the 
corresponding smallest value of c is also indicated. 

We see also that the values of c for a fixed required accuracy increases with g, using 
fixed a, and decreases with a for a fixed g. As expected, in any situation, c increases 
as e decreases. Regarding the values of Tg^c, we see that they form a monotonic series, 
increasing with respect to both g and c, approaching T from the lower side. In order to 
extrapolate the finite series results for a limit value T and z/, we have used the Vanden 
Broeck and Schwartz (VBS) extrapolation procedure used both in ([7]) and We 
observe that extrapolated values for z/ = limc,^oo ^g, which were estimated on the for Tg^c, 
get better for small e is reduced and c increases. 

In Table 5, we show the VBS extrapolated values of u for the same values of a and a 
fixed value of e that is small enough. Here we are lead to the most interesting result of this 
paper. For the purpose of comparison, we use the known exact values for u = l/(a — 1) 
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1/ 


y 


y 


1.2 


3.803537 


7.0 


5.0 


1.4 


2.437089 


2.7 


2.5 


1.6 


2.109820 


2.15 




1.8 


2.246300 


2.22 





Table 5: Comparative table of VBS extrapolated values of T and for distinct values 
of 1 < a < 2, using e = 10"'', obtained from the data in Table 1. We also include the 
corresponding value calculated by Glumac and Uzelac [3, z/^^, and the exact value of 
Vex for the mean field values of a. 

when 1 < a < 1.5, when the values obtained from the mean field analysis should prevail 
also for those from an exact solution. It is clear that our VBS extrapolation leads to a 
better agreement in comparison to the corresponding VBS extrapolated results in ([7]). 
Note that we use larger values of (A^ ~ 15 for all values of a) than those ones used 
in reference 0, which should be very important for the estimation of the the critical 
exponents. In this sense, the efficient TM method that we have applied is very useful to 
reduce the CPU time of the numerical calculations. 

The value of Vg for c = are indeed far away from the exact ones. This indicates that, 
although the value of T can be evaluated with good precision from the c = finite size 
data, the critical properties require indeed to probe the system in a close neighborhood of 
T. So, it is amazing to observe that our TM approach was able to be adapted to perform 
this most sensitive task. 

4 Conclusions 

In this work we presented results for the critical exponent v for the long range Ising 
chain based on the numerical evaluation of the eigenvalues of a 2 x 2 TM that condenses 
the information regarding the Boltzmann weights for all configurations of a finite size 
chain oi g + c + 1 spins including interaction among spins up to g sites apart. The 
adopted approach requires the minimum possible storage space and avoids the necessity 
of eigenvalue evaluation. The results were obtained with a double precision Fortran code 
implemented on a common desk computer. 

The comparison of our estimates with similar results reported by other authors shows 
that they are of the same quality or better than those obtained by a TM procedures that 
requires much larger matrices. In particular, we have used the known exact values of v 
in the range 1 < a < 1.5 to check the validity of our results. 

As compared to our previous implementation of our TM framework, we have shown 
that it is indeed reliable and that it could be successfully extended to much larger values 
of c, what was required in order to probe the close neighborhood of T. Despite the fact 
that analytical expression for the derivatives of the eigenvalues represents an important 
achievement of our approach and avoids the numerical differentiation of ^, we have ob- 
served that numerical overfiows in our Fortran code did not allow pushing the series of Vg 
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to the same larger values of g = 24 that were reached in our first evaluation of the critical 
temperature when c = 0. Efi^orts to sidestep this efi^ect and to include the evaluation of 
the magnetization exponent (3 are currently being undertaken to present a more complete 
analysis of this model. 
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